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Within the framework of the variational approach the ground state is studied in a gas of Fermi 
atoms near the Feshbach resonance at negative scattering length. The structure of the originating 
superfiuid state is formed by two coherently bound subsystems. One subsystem is that of quasi 
molecules in the closed channel and the other is a system of pairs of atoms in the open channel. 
The set of equations derived allows us to describe the properties of the ground state at an arbitrary 
magnitude of the parameters. In particular, it allows one to find a gap in the spectrum of single- 
particle Fermi excitations and sound velocity characterizing a branch of collective Bose excitations. 

PACS numbers: 03.75.Kk, 74.20.Mn, 67.60.-g 



Recently, the trend of investigations, associated with 
ultracold gases of Fermi atoms and Bose molecules gen- 
erated on their base, has been established and actively 
developed. The key point here is a use of the Feshbach 
resonance for the pair interaction of low-energy Fermi 
particles. The dependence of the resonance position on 
steady magnetic field B has given a unique opportunity 
for continuous variation of the effective interparticle in- 
teraction and for the change its sign (sign of scattering 
length a) with variation of magnetic field 

Several stages of studies in that field elapsed for a rel- 
atively short time interval. At the first stage there is 
demonstrated a production of ultracold molecules 40 K 2 
@ and 6 Li 2 0, 0, @ in the vicinity of the Feshbach 
resonance at magnitudes B corresponding to the posi- 
tive scattering length and, thus, existence of the gen- 
uine bound state near the dissociation threshold. In both 
cases the resonance occurs in the channel of S'-scattering. 
For fermions polarized over electron spin, this supposes 
the interaction of non-identical atoms in different hyper- 
fine states. 

It is found that the molecules produced live relatively 
long in such excited state. This supposes a possibility 
for the quasi adiabatic transition with variation of field 
B from the state with a < in which genuine molecules 
do not form into the state with a > and vice versa. 
The nature for suppression of inelastic processes during 
collisions with the transition of molecules to deeper levels 
is explained consistently in [||. 

The conservation of molecules at one vibrational level 
facilitates the conditions for Bose-Einstein condensation 
(BEC) in originating gas of composite bosons. In the ex- 
plicit form the condensation of excited molecules is dis- 
covered in a gas of atoms 40 K as well as in 6 Li 

The greatest interest in this field is connected with 
the study of Fermi atoms near the Feshbach resonance 
at magnetic fields corresponding to the scattering length 
a < 0. 

A possibility to choose the parameters, converting rar- 
efied gas into the system of strongly interacting parti- 
cles with the simultaneous lack of real molecules, makes 



an analysis of many-particle superfiuid state, which ap- 
pears in these conditions, especially interesting. For the 
first time, the question on the use of the Feshbach res- 
onance for achieving superfluidity in a ultracold gas of 
Fermi atoms is discussed in [l^. The problem 

of reconstructing the superfiuid state in the course of 
the transition from the weak attraction to the st rong 
one is considered earlier in the known papers 
The interaction had a single channel character and the 
spin configuration remained unchanged. The specific fea- 
ture of the Feshbach resonance lies in the two-channel 
character of the interaction. Particles in the continu- 
ous spectrum interact via the formation of the interme- 
diate quasi-molecule state in the other spin configura- 
tion. A snapshot fixes a two-component character of the 
system with pairs of particles in the quasi-molecular or 
unbound states. A dynamic coherent exchange realizes 
between two subsystems, providing permanency in equi- 
librium at the amplitude level of the population for the 
quasi- molecular state (decaying in vacuum). The coher- 
ent superposition of the quasi-molecular states and pairs 
of atoms in the continuous spectrum provides the forma- 
tion of the superfiuid state with the unified condensate. 

For the first time, the picture of such state is revealed 
experimentally at first in gas 40 K 15] and then in gas 
6 Li ^(|. The idea of experiment is in the adiabatic pro- 
jection of the a < state onto the genuine molecular 
state with a > due to sufficiently rapid variation of 
field B, preventing from the formation of the condensate 
in the course of the transition. Special measurements 
have confirmed the latter assumption, see, in particular, 
recent paper ^tJ- The experimental results in this field 
have stimulated the appearance of a series of theoreti- 
cal works El,[lii,|2(i,|22,[22 with the attempt of the 
microscopic description of the strongly interacting Fermi 
atoms, employing, however, one or another variants of 
the perturbation theory. 

The present work is devoted to finding the ground state 
of the interacting Fermi atom system near the Feshbach 
resonance at a < within the framework of variational 
approach. As a result, there appears an opportunity for 
unified description of the ground state at an arbitrary 
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relation of parameters. 

A choice of the macroscopic variational function 
reflects the fact that the ground state of the system in 
the vicinity of the Feshbach resonance is a unified state 
formed in essence by two coherently bound subsystems, 
namely, subsystem of Fermi atoms in the BCS-like state 
and subsystem of Bose quasi-molecules in the conden- 
sate state. In addition, ^f contains the states of quasi- 
molecules with nonzero momentum of the gravity cen- 
ter k ^ and those of pairs of atoms with nonzero 
total momentum. This allows from the very beginning 
to involve abovecondensate excited states (temperature 
T = 0), which play a noticeable role during the forma- 
tion of the ground state, into the variational procedure. 
The system of equations obtained gives an opportunity, 
in particular, to find a gap in the spectrum of single par- 
ticle excitations and velocity of sound characterizing a 
branch of collective Bose excitations and, thus, to deter- 
mine the critical velocity. 

In the present work we give a solution of this set of 
equations in explicit form for the limiting cases which 
allow an analytical solution. We restrict ourselves with 
an assumption of the smallness of the resonance coupling 
constant. The numerical solution in the general case will 
be given elsewhere. 

To make analysis more transparent, we consider the 
uniform system. The exchange between the subsystems 
of Fermi atoms and quasi-molecules with the conserva- 
tion of the total number of atoms determines the appear- 
ance of a single chemical potential fj,. The generalized 
Hamiltonian of the system can be written as 

H' = H-nN = ^^a+^ + ^tkHh + V + Va + Va- 

P,c k 

(1) 

Here a pa is the annihilation operator of Fermi atom in 
the open channel, a — 1,2 identifies the hyperfine state of 
an atom in the binary mixture and bk is the annihilation 
operator of quasi-molecules. 



interaction between atoms in the open channel in Q 



o 2^ a pi,i a -pi,2 a -P2,2a P2i i, 



VI = Ui 

pi,p 2 ,q^0 



(3) 



Si,l^g-pi,2 a g-Z>2,2%> 2 ,l- (4) 



The Hamiltonian in the form in fact, coincides with 
that widely used for the analysis of g ases under Feshbach 
resonance, see, eg., O , > > HIT H3l The sin S le dis- 
tinction is in the involvement of interaction which 
role can be noticeable in the general case. Strictly speak- 
ing, in Hamiltonian Q the terms with the nonresonance 
interaction of atoms with quasi-molecules and between 
quasi-molecules are omitted. However, both these inter- 
actions are produced by the resonance interaction J2J). 

The many-particle variational function is taken in the 
form 

*0 = *M*M*A |0> (5) 

The field operator of atoms has a usual representation 
= K + Vp^t^-pa] 



In the coherent representation the field operator of quasi- 
molecules in the state k = has the form 



e~2 M ° exp 



It is assumed that the number of such quasi-molecules 
M is large. 

A special feature of the structure for the many- 
particle variational function © is an involvement of the 
states connected with virtual production of excited quasi- 
molecules. The corresponding operator is represented as 

*M = II + Gk X! ^fc «fc/2-p,2afc/2+p,l] 
fc#0 p 



i p = ftV/2m - fi, 0c = Co + h 2 k 2 /Am, Co = £o ~ 2/i, 

where £o is a position of the resonance level in the closed 
channel, m is the atom mass (we neglect a difference be- 
tween hyperfine energies in states 1 and 2). The operator 



Note that the order of operators in (JSJ is essential. Op- 
erator acts on the the state resulted from the action 
of operator on the vacuum. The normalization of 
function 10 results in two relations 



1, F* + w k G 2 k = l 



(G) 



V - g^2(i>t 0-k/2-p^k/2+p,l + afe/2+p,l"fc/ 
p,k 



2-p,2 



b k ) (2) 



where 



w k 



E2 2 
V p+k/2 V p-k/2- 



(7) 



corresponds to the exchange resonance coupling between 
particles in the open and closed channels. Hereafter, we 
suppose that volume = 1. 

In the explicit form the coupling between pairs with 
zero total momentum is selected from the nonresonance 



We assume that coefficients u p , v p , i 7 ^, Gk are real. 

The variational procedure reduces to finding the en- 
ergy 



E = (*o H' * 
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and varying it over functions v p , F k , Mq with account of 
relations 10. Taking the structure of variational function 
© and the form of Hamiltonian JJJ into account, we find 



E = 2 E P C P «p + CoMq + 2VMo E P 
+^o(E p u p v p) +Uo(E p v p 



(8) 



For unit volume which we took, the magnitude Mq 
equals, in fact, to the density of quasi-moleculcs with 
k = 0. 

For the successive derivation of (jHJ, functions v p are 
multiplied by a factor exp(— z„) and u p v p by factor 
exp(— 2z p ) where z p = J2 k G\v p v p _ k . The direct calcu- 
lations show that magnitude z p proves to be small com- 
pared with unity, at least, in the region of parameters 
interesting for us. The involvement of these factors prac- 
tically has no influence on the results of the variational 
procedure. We put them equal to 1. 

Varying over v p , we find 



SE 
6v 



I = A^ p v p + (u p - ^) [2VMo.g + 2U Q £ p , u p > v p ,} 
E M o l^iCkGi + 2gG k F k - gw k %) = 



Varying over G k and taking © into account, one has 
SE 2w k 



(9) 



[CkG k F k + g- 2gw k G 2 k ] = 



(10) 



SG k F k 

Using this equation and the relation resulted from (JJJ 
dw k 2 

we rewrite @ in the form 
2£' p u p v p + (1 - 2v 2 p ) [y/Mo~g + U J2 W] = °- ( n ) 



Here 
6p = 



dp, dp = J] v 2 p _ k gG k F h + U J2 v l>- ( 12 ) 



Finally, varying over Mo, we arrive at 

£° + 9—rfr^u p v p = 0. 
-'AZn — 



(13) 



The extremum determined by conditions i|10|) . (|ll|l . and 
(I13|l corresponds to the energy minimum provided that 

gu p v p < 0, gG k F k < 

Hence, in particular, it follows that the first term in l|12|) 
and the second in (|13fl are negative. For definiteness, 
henceforth we put g > 0. 



Let us introduce the notation 

A = g\ZM~o~ +U '^2 u p v p 
p 

Solving equation (|11|) with @, we find the expressions 
similar to the BCS theory 



2 Up Vp 



A_ p 



Then the equation for the gap A reads 

p ^ 

At the same time equation l|13f) goes over into 



e Q -2n = 



- 1 V 



EV, 



(14) 



(15) 



(16) 



It is interesting that the role of virtual excited quasi- 
molecules in (|15|l and (|16fl is reduced to the renormaliza- 
tion of the atom spectrum, see, l|12[> and (|14|l . 

Solving equation (|10fl with JHJ allows us directly to find 



Gl 




(17) 



Equations (|15fl and l|l(jD contains formally an integral di- 
verging at the upper limit. This divergency is related 
with the character of introducing parameters sq, g and 
Uq into Hamiltonian JJJ . In the general case this requires 
the introduction an upper cutoff for the momentum re- 
gion region in which the interaction realizes effectively. 
However, in some important cases such cutoff arises in 
the natural way. 

Let us consider the most interesting case when the res- 
onance interaction Uq = is determining. Then from Eq. 
(11511 we have 



A = gy/M 



(18) 



The magnitude Eq in Hamiltonian (JTJ corresponds to the 
position of the resonance level. At g ^ the level renor- 
malizes even in the limit of density tending to zero, i.e., 
in vacuum. In this limit this can straightforwardly be 
traced by calculating the scattering amplitude for a pair 
of atoms with zero total momentum which is determined 
by interaction J3J) . Summing a series of the perturbation 
theory at E — > 0, we obtain, (cf. 



1(E) 



hs/y/m 



,3/2 



E-e + y 2 EpEp 1 +isy/E' 



9 



It is seen that the resonance level energy at g ^ renor- 
malizes and takes the value 

£a=e --g — 
z p £ p 
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In the many-particle problem it is natural to treat mag- 
nitude e as a real position of the resonance level. Sub- 
tracting \g 2 J2 p i from the both sides of iJTfljl and taking 

(118(1 into account, we have 



- „ 1 ,v^,l 1 , 



(19) 



Provided that the nonresonance interaction is involved, 
the situation becomes more complicated. However, treat- 
ing (|15|l and l|ltj|) jointly, we find the equation which again 
does not contain divergency 



1 

i -2ij,= -g 



U AA 



(20) 



If one employs the known relation between Uq and non- 
resonance scattering length a& fl beyond the Born approx- 
imation, see [2^], the gap equation H15|) can be trans- 
formed as 



A- 



1 . Anh 2 1 



(21) 



For g — 0, this equation reduces to the one familiar in 
the theory of superconductivity, see [13j. Assuming that 
the gas parameter is small for the purely nonresonance 
interaction, Uo in (|20() and (|12|) can be replaced with 
4irh 2 ab g /m. Solving equations (|2*U)l . J23J and (fTTfl to- 
gether with ©, {EJl ancl 031 , we find values A, M , G\ 
as a function of chemical potential \i for fixed (7, a^g, and 
£o- The values obtained allows us to determine the num- 
ber of atoms 



N A 



(*o| a+ a a pa |*o> = 2 £ t£ = £(1 - 

p,(T p p P 



and the number of virtually excited molecules 



(22) 



M' = (*o\b+b k |* ) =$>*<3 (23) 



The relation for the total number of atoms in the system 



N A + 2M + 2M' = N t 



(24) 



determines the equation to find chemical potential /i(N t ). 

As a result, the relative values N A /N t , M /N t , M'/N t 
can be found as functions of £o/2e t for the fixed dimen- 
sionless parameter 
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(25) 



Here Et is the Fermi energy corresponding to the total 
density N t of noninteracting atoms. It is worthwhile to 
note that the system of equations obtained holds only 



for the interval of variation So/2st in which condition 
M ^> 1 conserves. 

Let us make an important remark. In the strongly in- 
teracting many-particle system the superpositional char- 
acter of the ground state near the Feshbach resonance 
makes the definitions atoms and molecules relatively con- 
ditional. In fact, this is projections of the many-particle 
wave function onto the corresponding states. 

In the present paper we give a solution of the system 
of equations obtained in some limiting cases which allow 
analytical description. Let us first be g — > and Uo > 0. 
From Eq. (ED it follows that A = and from (TJJ it 
goes G fe = 0. 

According to (|20|l under these conditions the chemical 
potential fi — £q/2 for an arbitrary Uq =/= 0. The number 
of particles in the atom subsystem equals 



N A = 2£*(-£). 



with quantity £' p equal to 

in accordance with H12|l. Taking 124(1 into account and 
that M' o = 0, we find straightforwardly that in the 
interval 



< ie <£t + \u N t 



(26) 



gas of atoms in the normal state is in equilibrium with the 
condensate of quasi-molecules. At the left-hand bound- 
ary at Eq — > the number of atoms N A — > and all 
particles prove to be in the condensate of quasimolecules. 
As £o grows, the magnitude Mo falls continuously. So, 
at the right-hand boundary of interval (|26f) . N A = Nt, 
i.e., the purely atomic phase appears. For Uq < 0, the 
picture of behavior Mq and N A remains qualitatively 
same. However, the condensate of quasi-molecules van- 
ishes with growth of tildeeo before, approximately at 
ie — £t — I \ Uq \ Nt- In fact, this shift enhances to some 
extent owing to appearance of the gap A in the spectrum 
of atoms, which can be found from solving Eq. (|21l) at 

Let us turn to the solution of the system of equations 
obtained above when the resonance interaction is domi- 
nant. We omit the terms with Uo and restrict ourselves 
with values g which correspond to the small dimension- 
less parameter 7 <c 1. Then in accordance with (|18fl and 
one has A/e t = j^M /N t < 1. 

One can straightforwardly see that | d p \ ~ 7™ with 
m > 5/2 in 1(12(1 and the renormalization of the spectrum 
can be neglected. Let us start from the consideration of 
interval Eq ^> A. In this case from Eqs. ((19(1 we can find 



- 7 2 (^) 1/2 lnWA) 



(27) 
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Calculating Q22[l. we have for the number of atoms 



Na 
N t 



.3/2 



A 2 



Determining quantities Wk and G\ i|17[l , for the num- 
ber of abovecondensate quasi-molecules l|23|) in the limit 
£o > A we obtain 



~ I 3/2 /> 

iV t 2 7 U 



Employing relation (|24(l . we find the number of conden- 
sate quasi-molecules 



2M 



3/2 , / N 9/8 



(28) 



The magnitude (eo/^et),,, where Mo vanishes, one can 
estimate if the second term in (|27|) is omitted. 



2e t 



! 7 3/2 

3 r 



(29) 



In narrow interval of about 7 3 / 2 in the vicinity of this 
boundary the number of excited quasi-molecules M' 
proves to be larger than Mq. As e"o decreases, the num- 
ber of condensate particles M increases rapidly and 
M > M'. 

Let us find now the magnitude of sound velocity C in 
the system. For the equilibrium exchange of atoms be- 
tween both subsystems, the sound velocity is determined 
by the thermodynamic relation 



C 2 = 



m dNt ' 



(30) 



characteristic for the single-component systems, see, eg., 
|2^ . For g ^ 0, one has C 2 = 0. This result is a conse- 
quence of the transition of Fermi particles with compres- 
sion to the quasi-molecular subsystem and dp/ dp = 0. 

The similar effect takes place in an ideal Bose-gas at 
T <T C when a fraction of abovecondensate particles con- 
verts to the condensate with the pressure growth. 

When calculating the derivative dp/dN t , only the de- 
pendence A on Mo in the logarithm argument proves to 
be essential. (The involvement of the dependence of the 
second term in (|27|l on p, gives corrections ~ 7 2 ) 

djJ_ = _3_ 2 i /2 1 dM 
0N t ~ 16 7 [fi£t) M dN t 

Let us use the relation (|28|l omitting the last term on the 
right-hand side in the course of calculating dMo/dN t . 
In the derivative there appears a term proportional to 
A7jJ J^- , which plays a significant role near the boundary 
(|29|l within the interval in which Mo/N t becomes ~ 7 2 . 



Returning to definition (|30|l . we find a general expression 
for the sound velocity 



C 2 



CI 



1/2 



3mr<2 

277°* 



2 = J3_ 2ft (tV' 2 Nt_ 
' °* 16 7 m UJ M 



Far from the boundary (|29|l 

C 2 = C 2 

In the limit /i — > St and Mq — * wc find 

6 m 



(31) 
(32) 

(33) 



It is interesting that this limiting value is reached from 
the region where Mo ^ and the resonance interaction is 
absent. Thus the sound velocity as a function of io/2e t 
grows continuously to value (|3^|l typical for the collective 
oscillation of an ensemble of free Fermi particles. For 
7 2 ^ 1, the drastic growth C 2 to the value (|33|l occurs 
within narrow interval of about 7 2 . 

Consider now the region < eo < A. Equation l|19|l 
in this limit goes over into 



fi ~ -eo + 0.8g- 



(2m) 



3/2 



4tt 2 



-A 1 / 2 



(34) 



Accordingly, from (|22(l we find 



N t 1 V Nt 



3/4 



Determining the value M' from l|23(l . we have 
2M' 1 



N f 



:1 



21/8 



From the values obtained it follows directly that most 
of particles in this region are in the condensate of quasi- 
molecules 



2M 



1--7 
2 ' 



3/2 



1 



to 



21/8 



Employing relations (|30|l and (|34|) , we find for the sound 
velocity 



C 2 ~ 

3m 



(35) 



From l|34|l it is seen that the chemical potential tends to 
zero at 



£"0 



-7 5/2 £ t 



In order to obtain a complete picture, let us consider 
the region of negative e 0l assuming that A <C |e~o|- It 
is evident that real molecules dominate in this region. 
Admixture of atomic states is small. 



Na 
N t 



V 



2e t \ 1/2 M 



£0 



Nt 



(36) 
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The number of abovecondensate molecules is small as 
well. 

2M' 1 3 (2e t ^ 3/4 



N t 2 7 \s 



Correspondingly, the number of molecules in the conden- 
sate is 

~ l _ 1 2 (teA 1/2 _ 1 3 (2et} 3/4 
N t ~ 8 7 \e J 2 7 \e 

Equation ljl9|l in this case takes the form 



M = l: £ o + 0.55g 



2 m 3 / 2 |/x| 1/2 m 3 / 2 



+ 



2 ~ U ' "' ^ 2 32tt |e | 3/2 



9 4 M . 



Hence, involving l|25|) . we find for the sound velocity 



256 7 to V leo 



(37) 



The appearance of the high power of small parameter in 
(|37|l is due to the interaction between molecules, which 



in this case occurs via virtually created Fermi atoms. 
The number of the latter atoms is small, see l|36[) . Ac- 
cordingly, C 2 tends to zero at |eo| — ► oo. Comparing 
(|3T|) . and we see that the sound veloc- 

ity increases gradually with the growth £~o, reaching its 

maximum value jnnj. 

Superfluidity and critical velocity v c in the system con- 
sidered are determined simultaneously by the spectra of 
single-particle Fermi excitations and collective Bose ex- 
citations. Accordingly, v c takes the minimum value from 
two velocities: C and Tsxm{E p /p). Quantity E p is deter- 
mined by 1|14[1. For £o in the vicinity of 2e t , the critical 
velocity is governed by Fermi excitations and v c — ► 
together with M 0, A -> 0, see (JTSJ. Velocity 
min(Ep/p) increases monotonously as £o decreases. On 
the contrary, the sound velocity falls. The critical ve- 
locity therewith keeps to be governed by single-particle 
excitations up to £o/2e t ~ 1/3. For smaller values e"o/2et 
including the region of negative £o: the critical velocity 
is determined by the sound velocity. 

The work is supported by the Russian Foundation of 
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